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*5 ' Abstract 

We study the Polyakov loop and the correlator of two Polyakov loops at finite temperature in the 

d : 

weak-coupling regime. We calculate the Polyakov loop at order g 4 . The calculation of the correlator 

P^ | of two Polyakov loops is performed at distances shorter than the inverse of the temperature and 

>■ 

for electric screening masses larger than the Coulomb potential. In this regime, it is accurate up to 

order g . We also evaluate the Polyakov-loop correlator in an effective field theory framework that 

in 

r — | , takes advantage of the hierarchy of energy scales in the problem and makes explicit the bound-state 

O ' 

^^ | dynamics. In the effective field theory framework, we show that the Polyakov-loop correlator is 

at leading order in the multipole expansion the sum of a colour-singlet and a colour-octet quark- 
s' antiquark correlator, which are gauge invariant, and compute the corresponding colour-singlet and 
colour-octet free energies. 

PACS numbers: 12.38.-t,12.38.Bx,12.38.Mh 



I. INTRODUCTION 

The Polyakov loop and the correlator of two Polyakov loops are the order parameters of 
the deconfinement phase transition in SU(iV) gauge theories |l|, |2|. The phase transition is 
signaled by a non- vanishing expectation value of the Polyakov loop and a qualitative change 
in the large-distance behaviour of the correlation function (from confining to exponentially 
screened) |2j. In the deconfined phase, these quantities provide information about the electric 
screening and can be calculated at sufficiently high temperatures T in perturbation theory. 
For the correlation function of Polyakov loops, the validity of the perturbative expansion is 
limited to distances r smaller than the magnetic screening length r <C l/(g 2 T) [3|, |4|. 

From a phenomenological perspective, the Polyakov-loop correlator is interesting because 
it provides an insight into the in-medium modifications of the quark- ant iquark interac- 
tion. Indeed, in-medium modified heavy-quark potentials, inspired also by the behaviour of 
the Polyakov-loop correlator, have been used since long time in potential models (see e.g. 



Ref 



■B). 



However, although the spectral decomposition of the Polyakov-loop correlator is 



known, its relation with the heavy-quark potential is still a matter of debate and in need 
of a clarifying analysis [6[. The issue has become particularly relevant since recently an 



in-medium modified heavy-quark potential has been derived rigorously from QCD [7H11]. 
One of the aims of the paper is to discuss, in the weak-coupling regime, the relation between 
the Polyakov-loop correlator and these recent findings. 

The Polyakov-loop correlator is a gauge- invariant quantity, hence it is well suited for 



lattice calculations. In fact, the corre 



ator of two Polyakov loops has be calculated on the 



- |l5J as well as for full QCD [l6|, |l7| (for a review 



lattice for the pure gauge theory [12 

see Ref. |18|). Surprisingly, not much is known instead about the corre 

theory. The correlator is known at leading order (LO) since long time 

order, it was computed only for distances of the same order as the electric screening length 



ator in perturbation 
Il9| ; beyond leading 



in Ref. [20 1. 



The purpose of the paper is to evaluate the (connected) Polyakov-loop correlator up to 
order g 6 at short distances, rT <C 1. This corresponds to a next-to- next-to-leading order 
(NNLO) calculation, if we count the order g 4 as LO and the order g 5 as next-to-leading order 
(NLO). We also revisit the calculation of the expectation value of the Polyakov loop at order 
g 4 , which corresponds also to a NNLO calculation, if we count 1 as the leading-order result 



and g 3 as the NLO one. We will find a result that differs from the long-time accepted result 
of Gava and Jengo [2l| . Finally, we will add on the discussion about the relation between 
the Polyakov-loop correlator and the in-medium heavy-quark potential. 

The paper is organized as follows. In the next section, we discuss the gluon propagator 
in static gauge at one-loop level. Section UTTI contains the calculation of the Polyakov loop at 
NNLO, while in section HVl we calculate the Polyakov-loop correlator. In Sec.|Vj we rederive 
the Polyakov-loop correlator in an effective field theory language. There, we also define a 
singlet and an octet free energy that we compute. Finally, section IVT1 contains the conclusion 
and outlook. 

II. THE STATIC GAUGE AND THE SELF ENERGY 

The Polyakov loop and the Polyakov-loop correlator are gauge- invariant quantities. We 
may exploit the gauge freedom by choosing the most suitable gauge. A convenient gauge 



choice is the static gauge 22], defined as 1 

d A°(x) = 0. (1) 

The reason for using the static gauge is that in this gauge the Polyakov line has a very 
simple form 

L(x) = P exp (ig f^ drA°(x, r) J = exp (^^j , (2) 

where P stands for the path-ordering prescription. The spatial part of the gluon propagator 
reads 

A> n , k) = 1 [8 i3 + ^) (1 - M + ^ (s i3 - (1 - 0^) 5 n(h (3) 

where u n = 2nTn, n e Z, are the bosonic Matsubara frequencies and k 2 = bj\ + k 2 . 
Throughout the paper italic letters will refer to Euclidean four-vectors and bold letters to 
the spatial components. The parameter £ is a residual gauge-fixing parameter. We call 
non-static modes those propagating with nonzero Matsubara frequencies and conversely we 
employ the term static mode for the zero mode. The first term in the r.h.s. of Eq. ([3]), 
proportional to (1 — <5„o), is then the non-static part, whereas the second, proportional to 



1 We will work in Euclidean space-time and will label the Euclidean-time component. 



5 n o, is the static part. The temporal part of the gluon propagator reads 

Ax>K,k) = ^, (4) 

which is purely static. Note that the gauge-fixing parameter affects only the static part of 
the spatial gluon pr opagat or. The complete set of Feynman rules in this gauge has been 



discussed in Refs. 22h24|. Feynman rules are listed in appendix [A] together with our 
Feynman diagram conventions. We will adopt the static gauge in all the calculations of the 
paper, if not otherwise specified. 

A necessary ingredient for the calculation of the Polyakov-loop expectation value and the 
Polyakov-loop correlator at NNLO is the temporal component of the gluon self energy at 
LO. In the static gauge, due to the static nature of the temporal propagator in Eq. (@| only 
Iloo(k) = Iloo(0,k) enters. Furthermore, at LO static and non-static modes do not mix in 
Iloo(k), which can thus be conveniently split into 

LVk) = < s (k) + r4(k) + r4(k), (5) 

where the three terms correspond to the contribution of the non-static gluons, the static 
gluons and the fermion loops respectively. 




FIG. 1: Diagrams contributing to the non-static part of the gluon self-energy in the gluonic sector. 
Dashed lines are temporal gluons, curly lines are spatial non-static gluons. 

1. F4 s (k) 

In the gluonic sector, the non-static part of the self-energy receives contributions only 
from the two diagrams shown in Fig. [TJ Using the Feynman rules of appendix |Aj it 
can be written in terms of five dimensionally-regularized master sum integrals 

< s (k) = -2g 2 C A {^h - (d - l)h + I 2 + \h + ~/ 4 ) , (6) 



where Ca = N = 3 is the number of colours, d = 3 — 2e is the number of dimensions, 
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q = k — p, is a shorthand notation for the non-static, n 7^ 0, sum integral: 
/p 






2r 



(2vr) a 



(8) 



and /x is the scale in dimensional regular izat ion. The result fl6]) can be conveniently 
cast in a sum of a vacuum part, a matter part, a part made of the subtracted zero 
modes and a part that we may call singular, because it is singular for T — y 0; the 
singular part is a peculiar feature of the static gauge. We then have 
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(12) 
(13) 



where je is the Euler constant and 71b(&) = 1/ (e fc / T — l) is the Bose-Einstein 

distribution. We refer the reader to appendix [B] for details on the derivation of 

;hese equations. The vacuum part (fTU|) agrees with the static gauge computation in 

23| . Furthermore, the vacuum part and the matter part are identical to the k° — y 



imit of their Coulomb gauge counterparts, computed respectively in [25|, |26[ and 



27 



28]. n^ s (k) zero consists of the subtracted zero modes. In the e — y limit, it is 



T|k|/4; we have kept the order e corrections, because, in the Polyakov-loop correlator 
calculation of Sec. IIVt we will need to evaluate the Fourier transform of |k| 1_2e /|k| 4 , 
coming from a self-energy insertion in a temporal-gluon propagator, which is divergent. 



2. n£ (k) 

At leading order in the coupling, IlQ (k) may be written in terms of three dimen- 



sionally-regularized master sum integrals jz 



nS,(k) = 2g 2 n f 



-I + 2h + -h 



(14) 



where 
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(15) 



q = p + k and u)„ = (2n + l)7rT are the fermionic Matsubara frequencies and nj is the 
number of massless quarks contributing to the fermion loops. Since no fermionic Mat- 
subara frequency vanishes, fermions are purely non-static. The fermionic contribution 
can be cast into a sum of a vacuum and a matter part: IlQ (k) = IIq (k) var . + Ilnr, ( k)^ ^ . 
After the Matsubara frequencies summation, the matter part can be read from 29 1 



C(k) r 



2tt 
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|k| + 2|p| 



2|p| 



(16) 



where n-p(k) = 1/ [e k ^ T + l) is the Fermi-Dirac distribution. The vacuum part is 
given by 



n£)(k), 



2L.2 



2^k 



3(4tt) 



i n f 



7e + ln(47r) 



, k 2 5 

In— + 



(17) 



e ^ 3 

3. n^(k) + n*(k) 

Let us now consider the sum Il^ s (k) + IiQ (k). The divergences in the vacuum parts 
(TTOj) and (TT7|) are of ultraviolet origin and are accounted for by the charge renormal- 
ization. In the MS scheme, the renormalized sum of vacuum parts reads 



21,2 



noo (k)vac + n 00 (k), 



g*k 



(Air) 



' , fi 2 31 10 



(18) 



where (3 = llC A /3 - 2n f /3. 

Simple analytical expressions can be obtained for the renormalized sum Il^ s (k) + 
IlQ (k) in the two limiting cases |k| <C T and |k| ^> T. In the former case, we have 



M s + r4)(|k|«r) 



2t^2 



9 2 T 



(c A + 



TV 

2 



g 2 k 2 [11 / , 4vrT 2 1 n 

— C A -ln^— ^ + l + 2 7s 

3 V /^ 



(4tt 



2 / (4ttT) 2 



n f ( -ln v 2 ; - l + 27£ + 41n2 



+ , 2 k 2 (£) , (19) 

where the leading-order term is momentum independent and can be identified with 
the (square of the) Debye mass thd, 

™ 2 ^ 9 -it( n+ t)' (20) 

which provides, in the weak-coupling regime, the inverse of an electric screening length. 
We note that Eq. (fT9~l) presents a logarithm of the renormalization scale over the 
temperature rather than over the momentum: this happens because in the limit |k| -C 
T the matter part produces a term proportional to k 2 /3oln(T 2 /k 2 ) that combines 
with the logarithm in the renormalized vacuum part ffT8|) to cancel its momentum 
dependence. 

In the opposite limit |k| ^> T, we have 

W + C) (|k| » T) = < s (k) vac + C(k) vac + g*C A f-^ - J^) 

+g 2 C A M 4 M [1 + e(-l - lE + In(167r))] 
+9 2 T 2 n^\ . (21) 

We observe that, in this limit and at the considered order, fermions enter only through 



their contribution to the vacuum part. It shoulc 
— g 2 CAT 2 /18 term appears also in Coulomb gauge [9|. 



be also noted that, while the 
27J, the term proportional to |k| 3 



is instead a peculiar feature of the static gauge. The terms proportional to eT|k| 1_2,E , 
which appear in the second line, come from the subtracted zero modes and contribute 
only when plugged into divergent amplitudes. Details on the derivation of these ex- 
pressions can be found in appendix O 

4. n s (k) 

The diagrams contributing to the static part of the gluon self energy are shown in 
Fig. [2j They are not sensitive to the scale T, since, by definition, static gluon propa- 
gators are just made of zero modes, however they are to the scale m D . Hence, when 



'' \ 




FIG. 2: Diagrams contributing to the static part of the self-energy: the dashed lines are temporal- 
gluon propagators, the wavy lines are static spatial-gluon propagators. Loops made of two static 
spatial-gluon propagators and of ghosts vanish. 

evaluating the static contribution, it is important to keep in mind that, if the incoming 
momentum is of the order of the Debye mass, then insertions of gluon self-energies 
of the type of Eq. ( TT91) into the temporal-gluon propagator need to be resummed 
modifying the temporal-gluon propagator into 



»*(«»>*) = uMjT- (22) 



k 2 + m 2 



The static part of the gluon self energy with resummed propagators reads, for all 
values of the gauge- fixing parameter £, 

^s /,n i„ m i f f d d p ( 1 d-2 

I4(k) = g 2 C A T^ / -£- -^—T + 



., (2n) d \p 2 + m 2 D p 2 
, 2(m 2 D -k 2 ) 2 2 p 2 + 2p ■ k 

+ 2? 2 i 2T + (^ -1 )( k + m D) 4/ 2 , H 

p 2 (q 2 + m|)) p 4 (q 2 + m|)) 

where q = k + p. The result agrees with Ref . J3|, [30| . Note that Eq. ( 125|) applies for all 
gauges sharing the same static propagator, among which the static and the covariant 
gauges. The expression is finite in three dimensions and reads 



n s fkl _ 9 2 C A T 



m o — k 2 k . . 

! — arctan thd + (c — 1)tid 

k rra^i 



(24) 



Finally for the static part |k| 3> T implies |k| ^> mp and 



ngo(|k| » m D ) = -/CU <! T|k| , ^ ' [1 + e(-7B + ln(167r))] + 0(m D T) \ , (25) 



4 
where again we have kept up to order e terms proportional to T|k| 1_2e . 

5. n 00 (k) 

iloo(k) is obtained by summing (|ID1>, (DD, (H2D, (II3D, <HSJ», (HZH and (|23D (or (J23D)- In 

8 



particular, the asymptotic expression for the gluon polarization at high momenta is 

g 2 k 2 ( /i 2 31 10 \ 2 / T 2 |k| 3 \ 



n o(|k|»T) = -^(/5olng + ^-^/) + ^(-^ 



192T 



VCa ' ' " + [9 2 TJ,9 2 mnT ) . (26) 



4 V k 2 

Note that the term proportional to T|k| 1_2e e° in Eq. ( 125]) has canceled against the 
term proportional to T|k| 1_2<E e° in Eq. ( I2T1) . 

III. THE POLYAKOV LOOP 

The quantity we are interested in computing is the trace of the Polyakov line L = Lr 
in a representation R of dimension d(R), where R is either the fundamental representation 
(R = F, d(F) = N) or the adjoint representation (R = A, d(A) = N 2 - 1): 

( iR ) = (f riR ), 4 = -^. (27) 

The brackets stand for the average in a thermal ensemble at a temperature T. Expanding 
the Polyakov line in the static gauge up to order g 4 yields 

/ r \ - 1 g2 ^ rA o) / ^ A D . S 4 (TrA 4 ) 

(L^) - 1 - — -j^ i— t?j + — TA +... . [2X) 



In computing Eq. ( 1281) perturbatively, each diagram can receive contributions from both 
scales T and mp, for which we assume a weak-coupling hierarchy: 2 

T > m D . (29) 

In the weak-coupling regime, the calculation of (Lr) may be organized in an expansion in 
the coupling g\ our aim is to compute (L R ) up to order g 4 . Sometimes, we will find it useful 
to keep tud/T as a separate expansion parameter with respect to g, in order to identify 
more easily the origin of the various terms. We will call the g 3 term the NLO correction to 
the Polyakov loop and the g A term the NNLO correction. We will also identify the source of 
some higher-order corrections of order g 5 and g A x (m^/T) 2 that will play a role in Sec. |VJ 

2 When discussing energy scales, we will consider T and multiple of ttT to be parametrically of the same 
order. 



V ' V / \ / 

V V V 



FIG. 3: Diagrams contributing to the perturbative expansion of g 2 (TvAQ). The dashed line is a 
temporal-gluon propagator, the dot is the point x where the loop originates. The blob stands for 
the gluon self energy. 

A. The order g 3 contribution 

Let us start examining g 2 {TiA^). Diagrams contributing to g 2 (TiA^) are shown in Fig. [3J 
Summing up all these diagrams, g 2 (TrA 2 ] ) can be written as 

6{Lr) ~ ~%~^~ ~ ~^F^ J (27r)"k2 + IIoo(k)' (30) 

where Cr is the quadratic Casimir operator of the representation R {Ca = N, Cp = (N 2 — 
1)/(2N)). We observe that the integral receives contributions from the scales T and mo- 
We set out to separate the contributions from these two scales assuming the hierarchy (1251) . 

1. Modes at the scale T 

We evaluate the integral (1301) for |k| ~ T ^> mo- In this momentum region, n o(|k| ~ 
T ~^> mo) "C k 2 and we may expand the gluon propagator in n o- The LO term yields 
a scaleless integral 

„2 r Adh. i 

(31) 



(32) 
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{2n) d 
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whereas the following term gives 
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d d k n 

(2yr) d 


oo(|k| 


k 4 


m D ) 


This term 


is of order g 4 . 













2. Modes at the scale itld 

We evaluate now the contribution from the scale m/> We recall from Eqs. (fTT?|) and 
(1241) that, for |k| <C T, Iloo(k) = m 2 D {\ + O(g)). We then rewrite the propagator in 

10 



Eq. flU asl/(k 2 + n 00 (|k| < T)) = l/(k 2 + m 2 D + (n 00 (|k| < T) -m 2 D )) and expand 

in IIoo(|k| <T)- m|,. The LO term yields 

e /r x 9 2 Cr 2e f d d k 1 C R a s m D 

S{L R ) LOmD = -^rV J {27r)dk2 + m 2 D = "a— jr. ( 33 ) 

whereas the following one gives 

, /r v 9 2 Cr 2e f d d k n 00 (|k|~m 1) «r)-m 2 ) 

WnU)bd " "^ il^p W+™& 2 ' ( } 

which is at least of order g A . 



Up to order g 3 , we then have 



<!«.> = l + %^? + 0(9 4 ). (35) 



/ > 
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/ 



A 



a) b) 

FIG. 4: Diagram a) is the leading-order contribution to (TVAq): it vanishes because of the three- 
gluon vertex involving only temporal gluons. Diagram b) is the LO term of g^^TicA 4 ,): it vanishes 
because scaleless. 

The LO contribution to the cubic term g 3 (TrAQ) is shown in Fig. H]a). It vanishes due 
to the structure of the three-gluon vertex. This is just a LO manifestation of the charge- 
conjugation symmetry; in fact, due to this symmetry, g 3 (TiA 3 ) ) vanishes to all orders. The 
quartic term g 4 (TrA' i l ) gets its LO contribution from the diagram shown in Fig. H]b), which 
vanishes because scaleless. At higher order, a comparison with the analysis we have just 
performed for (TtAq) makes it clear that g A {TrA 4 l ) starts to contribute at order g 4 x (m[)/T) 2 , 
which is again beyond the accuracy of this analysis. We can therefore identify as the only 
contributions to the Polyakov loop at order g 4 the ones of Eqs. (1321) and (1341) . In Sec. IIIIB1 
we will compute these contributions and, in Sec. IIIIDl we will analyze some sub-leading 
terms. 

11 



B. The order g 4 contribution 



We now set out to compute Eqs. (|32j) and (134"]) . Following the discussion in Sec. 011 
we separate the non-static from the static modes in n o(k). We then have four sources of 
contributions: non-static modes at the scale T, non-static modes at the scale rriu, static 
modes at the scale T and static modes at the scale rriD- 



1. Non-static modes at the scale T 
The non-static contribution to Eq. 
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TT n J (2n) d k 4 

where Il^ s (|k| ~ T) is the full non-static contribution as defined in Eq. ([9]) and 
similarly TTq ( | k | ~ T) is the full fermionic contribution as defined in Eqs. (TT6J1 and 
(IT7|) . We can rewrite Eq. (136]) as 
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where we have defined the two-loop master sum-integrals Jj and Ji as 
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These integrals are evaluated in appendix |D] and their sum yields 
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(39) 



The divergence stems from the J2 integral and is expected to cancel against an opposite 
divergence coming from the scale mo- 



2. Non-static modes at the scale tud 

The non-static contribution to Eq. ( 134]) reads 

d d k nNS(|k|~m D )+ng,(|k|~m D ) 



<*<£*> 
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For |k| much smaller than the temperature, Eq. (1T91) applies and thus Il^ s (k) + 
IToo(k) = m 2 D + (9(g 2 k 2 ). Therefore, the contribution of Eq. ( 140]) is of order g A x 
(mo/T) ~ g 5 . More explicitly, plugging Eq. f lT9|) into Eq. ( 140]) gives 



*<£*> 



NS,m D 



3g 4 C R m D 
4(4tt) 3 T 



/So 1b (jJ^) 2 + 2/3 o7£ + yC A - ~n/ (4 In 2 - 1 



• (41) 
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Although a term of order g 5 is beyond our accuracy, the contribution (1411) is of interest 
because it fixes the renormalization scale of g 3 in the LO term fl35l) (a s rriD/T ~ g 3 ) 
to [i — 4nT. 

3. Static modes at the scale T 

The static contribution at the scale T to Eq. (I3"2"j) reads 

It vanishes because IlQ (|k| ~ T ^> mc) ~ 5 ,2 ^|k| (see Eq. (125]) ) and thus the resulting 
integration over k is scaleless. 

4. Static modes at the scale mo 

The static contribution to Eq. (1541) is 



« Bfl = ^V* / ,0^,1.2.1.2^ ^ 



^Cfc 2e r d d k n s (|k |) 

2T ^ 7 (2vr) d (k 2 + m D 
where IIq ( | k | ) is the full static contribution of Eq. ( 1231) . The computation is carried 
out in detail in appendix [EJ the result reads 

KL R)smD = g^i (-1 - ,n £^ - I + 1E - 1„(4,)) . (44) 

The divergence cancels against the one of Eq. (1391) coming from non-static modes 
at the scale T. 3 Note that the gauge-dependent part of Eq. ( |23l) gives a vanishing 
integral, thus yielding the expected gauge-independent result. 

Summing all contributions (static and the non-static) from the scales T and tud up to 
order g A thus gives 

+ 0(g 5 ). (45) 



C R a s m D C R a 2 s 


Ca 


6CK 1 


2 T 2 




V T 2 2 



(LR) = 1 + -7T-£- + -^ C M'n^ + - -^2 



3 Both divergences in Eqs. (|39|) and (j44|) are of ultraviolet origin. This seems to contradict the expectation 
according to which infrared divergences from higher scales should cancel against ultraviolet divergences 
from lower scales. The contradiction is only apparent. The static modes at the scale T develop both an 
ultraviolet and an infrared divergence that cancel against each other if regularized by the same cut off in 
dimensional regularization as assumed in Eq. (|42j) . In general, however, the ultraviolet divergence of the 
static modes at the scale T cancels against the ultraviolet divergence of the non-static modes, such that 
the sum of static and non-static modes at the scale T ends up having only a residual infrared divergence. 
It is precisely this infrared divergence coming from the scale T, formally identical to the divergence in Eq. 
(|55)l . that cancels against the ultraviolet divergence in (|H|| coming from the scale m_o. 
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C. Comparison with the literature 

At order g 4 , the Polyakov loop was first calculated in the pure gauge case (n/ = 0) and 
in Feynman gauge, by Gava and Jengo (GJ) 21], who find 



C R a s m D , C R C A a 2 s f, m 2 D 3\ -. 



(^.)G J = l + ^^ + ^y^( v ln^-21n2 + ^+0(, 5 ). (46) 

Their result disagrees with ours, given in Eq. ( 1431) . 

The disagreement may be traced back to an incorrect treatment of the static modes at 
the scale mn in [21] . In Feynman gauge, at order g 4 , three terms contribute to the Polyakov 
loop: the non-static gluon self energy, whose dominant contribution comes from the scale T, 
the static gluon self energy, getting contributions from the scale m^ only, and a third term 
coming from the fourth-order expansion of the Polyakov line. The computation of Gava 
and Jengo correctly reproduces the first and the third term. We show this with some detail 
in appendix |F] However, in the evaluation of the static gluon self energy, the Debye mass 
is not resummed in the temporal gluons, leading to an inconsistent treatment of the scale 
mi). 4 Indeed, they have 

tt«? /, n 9^ nr, 9, f d d p fd-l 2k 2 \ 

l4(k) GJ = iC A Tf- J^—-—y (47) 

which is the static self energy in Feynman gauge but without resumming the Debye mass 
in the internal propagators. If, instead, the Debye mass is resummed, the expression of 
the static self energy changes to Eq. (|23l) with £ = 1. In this case, the calculation of the 
Polyakov loop in Feynman gauge leads to exactly the same result as in Eq. (1451) . 

While the last part of this paper was being completed, Burnier, Laine and Vepsalainen 



31] published a perturbative analysis of the singlet quark-antiquark free energy. In the first 
part of their work, they consider also the Polyakov loop at order g 4 within a dimensionally 
reduced effective field theory framework in a covariant or Coulomb gauge. Our result (14"5]) 
agrees with theirs. 



4 



In 



21 



some contributions coming from the resummation of the Debye mass seem to have been included 



m 



SW(0). 
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D. Higher-order contributions 

In Sec. IIIIBt we obtained in Eq. ( 14"T|) a term that is of order g A x [mr,/T) ~ g 5 . Other 
contributions of order g 5 can only come from (TtAq). Hence, they are encoded in the two- 
loop expression of the gluon self energy. 

At order g 6 , we can expect other contributions from the two-loop self energy and contri- 
butions coming from the diagram in Fig. H]b). We explicitly calculate these last ones due 
to their relevance for Sec. |V] The computation is carried out by evaluating the colour trace 
of the diagram in the representation R, whereas the loop integrations are easily obtained by 
comparison with Eq. (13"5]) . Thus we obtain 



S(L R ) 



{^- c -¥-)im 2 - w 



The colour structure of this quartic term is not linear in Cr, a fact that will play a role in 
Sec. |V] We recall here that the linear dependence of ln(Lpj,) on the Casimir operator Cr is 
called Casimir scaling of the Polyakov loop. Equation (HHj) provides the leading perturbative 
correction that breaks the Casimir scaling. It is a tiny correction of order g 6 , which may 
explain, at least in the weak-coupling regime, the approximate Casimir scaling observed in 



lattice calculations 32]. 



IV. THE POLYAKOV-LOOP CORRELATOR AT ORDER g Q FOR rT < 1 

The spatial correlator of Polyakov loops in the fundamental representation is defined 



as 

(Tr4,(0)TrL F (r)>. (49) 



Following the notation of [20|], we define Cpl(V, T) as the connected part of the correlator 
C PL (r, T) = (Tr4(0)frL F (r)) c = (Tr4(0)fr J L F (r)) - (L F ) 2 . (50) 

Expanding Eq. (1501 up to order g & yields 5 

r (r T ^ 9' (Tr^(0)fr^(r)) c g* (fr^(0)fr^(r)) c 
C PL [r,l) - — — — + 



(2!)2 T 4 ^ 3! )2 T e 



5 We adopt a slightly different definition of Cpl {r, T) with respect to [20( , in that we consider the zeroth- 
order term in the perturbative expansion, i.e 1, as part of (Lp) 2 rather than of Cpl- 
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2g« (fr^(0)fr^(r)) c . . 



G(g»). (51) 



2! 4! T 6 

Since the generators of SU(iV) are traceless, the first term in the expansion, which is 
g 2 (TrA (0)TrA (r)) c , vanishes and thus the correlator starts in perturbation theory with 
a two-gluon exchange term. Terms with an odd number of gauge fields have been omitted 
from Eq. (|5"T1) since they vanish for charge-conjugation symmetry. 

We will perform a complete calculation of the Polyakov-loop correlator for distances 
rT <C 1. This situation corresponds to temperatures lower than the inverse distance of the 
quark- ant iquark pair, hence it is the right one to make contact with known zero-temperature 
results. We assume the following hierarchy: 

1 Q 2 

- > T > m D > — . (52) 

typical scales appearing in any perturbative static quark- 



The scales 1/r and g 2 /r are the 

antiquark correlator calculation 33|, |34j . The scales T and mo are associated to the ther- 
modynamics of the system. We assume that they are smaller than 1/r, because we are 
interested in short distances. We assume that they are larger than g 2 /r, because we would 
like to study a situation where both thermodynamical scales affect the quark-antiquark po- 



[9 



tential [9|. In the weak-coupling regime, as discussed above, T 3> m D , where mo is given 
by Eq. ( |20l) . Equation ( 1521) amounts to having two largely unrelated small parameters, g 
and rT, the hierarchy only requiring rT ^> g. Differently from the Polyakov-loop calcu- 
lation where we had only g, the perturbative expansion of the Polyakov-loop correlator is, 
therefore, organized as a double expansion in g and rT. We will stop the expansion for the 
Polyakov-loop correlator at order g e (rT)°, meaning that, given a term of order g k (rT) n , we 
will display it only if k < 6, for any (positive or negative) n, or if k = 6, for n < 0; we will not 
display it elsewhere. We should note here that, as in any double expansion whose expansion 
parameters are unrelated, undisplayed terms may, under some circumstances, turn out to 
be numerically as large as or larger than some of the displayed ones. 6 

In 20| , Nadkarni computed the Polyakov-loop correlator up to order g 6 using resummed 
temporal-gluon propagators throughout the computation, which amounts to calculating the 
Polyakov-loop correlator for distances rrri£, ~ 1. Our calculation will differ from Nadkarni's 



6 A posteriori (see the final result in Eq. <j67j) ) , this may be avoided, in our case, by further requiring that 
rT » V5- 

1(3 



one in that we adopt the different hierarchy (j52|) . Nevertheless, some of our results can be 
obtained by expanding Nadkarni's result for rmo <C 1; we refer to Sec. IIV HI for a detailed 
comparison between the two results. 

The calculation of the different contributions to Eq. (ISTj) will proceed similarly to the 
calculation of the Polyakov loop performed in the previous section. We will consider the 
different Feynman diagrams contributing to each of the terms in (}5Tj) . separate the contribu- 
tions from the different energy scales and, in case, distinguish between static and non-static 
modes. 



V 

I 



V 

II 



V 

III 




FIG. 5: Diagrams contributing to (tr^(0)tr^(r)) c 



A. The leading-order contribution: diagram I 



We start by evaluating the four-field correlation function: its leading-order contribution 
is given by diagram I in Fig. EJ It does not vanish only for momenta of order 1/r, giving 

d d k e ~ ik - r ^ 2 



N 2 - 1 n A , 



8N 2 T 2 



/' 



(2iv) d k 2 



N 2 -l aj 

8N 2 ( rT )2- 



(53) 



B. The contribution from diagrams of type II 



As we go beyond leading order, the first class of diagrams that we consider are those with 
gluon self-energy insertions in one temporal-gluon line, whose first example is diagram II in 
Fig. They give 

1 1 



N 2 -1 g 4 1 
^p L (r,T) n = 2^^-— /, 



r\ d h 
2, / " K g-ik-r 



(27T) 



k 2 + n 00 (k) k 2 



(54) 
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where the factor 2 comes from the symmetric diagrams and Hoo is the sum of bosonic and 
fermionic contributions to the gluon self energy, as in the Polyakov-loop case. This diagram 
receives contributions from all scales and depends on the gauge parameter £. However it 
can be shown that the gauge dependence cancels with diagram IV [20], so, for simplicity, 
here we write our results in static Feynman gauge, £ = 1. 

1. Contribution from the scale 1/r 

We start by evaluating the contribution from the scale 1/r in the integral. If |k| ~ 
1/r > T, then we have 



5Cp L (r,T)ni/r 



N 2 -lg 4 1 
AN 2 T 2 ^ 



,2c 



/'" 



d d k _ jk . r n 00 (|k|>T) 



(27T) 



o 



rT 



(55) 



where Hoo(|k| 3> T) is given by Eq. (|26|) . The Fourier transform of the vacuum part 
corresponds to the one-loop static QCD potential and can be read from [35l. |36| . Using 
that in dimensional regularization the Fourier transform of l/|k| n becomes 37] 

d d k e" ikr 2- n 7r- d / 2 T (d/2 - n/2) 



we have 

SC PL (r,T) 1I1/r 



(27T) 



iV 2 -l a s 3 

r,, T2rT" 



„d—n 



r (n/2) 



(56) 



1 

2^ 



2/3 (ln(/ir) + Ts ) + — C^ - — n/ 



rT 7r(r 

9 



T) 2 )} + 0(^ 6 (rT) 2 ,^). 



(57) 



The term in the first line comes from the Fourier transform of the vacuum contribution, 
whereas the terms in the second line come respectively from the singular part, the (zero 
mode) order e term 7 and the T 2 term in Eq. (1261) . Higher-order corrections to Eq. 
(12"6"|) contribute at order g 6 (rT) 2 or g 7 . Higher order radiative corrections to the gluon 
self energy contribute at order g 8 . Note that the (a^/n)PQ ln(/xr) term in Eq. (|57"|) 
fixes the natural scale of a 2 in the LO term SCpi,(r, T)i to be 1/r. 

2. Contributions from the scales T and m^ 

We now consider the contributions from the thermal scales. For what concerns the 



7 The dimensionally-regularized Fourier transform of the order e term in Eq. (|26p yields a 1/e pole, 
eventually leading to a finite contribution. 



18 



temperature, |k| ~ T translates into r|k| <C 1 and nip <C |k|. Integrating out the 
temperature leads to the following contribution 



SCp L (r,T)n T 



N 2 -lg A 1 
AN 2 T 2 !^ 



// 



2c 



J^^O^.r) 



n 00 (|k|~T) 



x 



[i + oar 



(58) 



where we have implemented the condition r|k| Cl by expanding the Fourier exponent. 
Integrating out the Debye-mass scale leads to the following contribution 



SC PL (r,T) u 



N 2 -lg 4 



HID 



AN 2 T 2 Anr 



fJ 



2r 



d d k 
(2vr) d 

n 00 (|k| 



[l + 0((k-r 



k 2 + m 2 T 



D 



m 



D 



(k 2 + m 2 D ) 



2 \2 



+ o 



9' 



m 



D 



(59) 



The integrals to be evaluated are the same needed to evaluate Eqs. fl32|) . fl33l) and 
(!M|) . Thus, summing the T and m D contributions, we obtain 



5Cp L (r,T)ii T+ 



m D 



N 2 - 1 a 2 
AN 2 rf 
9 7 



m D 



T 



+ a s 



CL 



ln^ 



1 



T 2 + 2 



— n f In 2 



+C»(^.,^(rT) 



(60) 



The term of order g 5 /(rT) comes from the first term in (159]) . the terms of order g 6 /(rT) 
come from the non-static modes in fl58l) and from the static ones in the second term of 
Eq. (159]) . the appearance of the logarithm In m 2 D /T 2 signals the cancellation between 
divergences at the scale T and mo, the suppressed term g 1 jirT) comes from the non- 
static modes in the second term of Eq. (J5HJ) (see Eqs. (HO"]) and (HT]) for the analogous 
case in the Polyakov-loop calculation), whereas the suppressed term g 6 (rT) comes 
from the (k • r) 2 term in Eq. (158]) . 



C. The contribution from diagrams of type III 



Diagram III in Fig. [5] is the first example of the class of diagrams with gluon self-energy 
insertions in both temporal-gluon lines. They may be evaluated from the diagrams of type 
II: 



6C PL (r,T)ui 



N 2 -l g 4 
8N 2 f 2 
8N 2 T 2 

N 2 -l~g~* 



[J 



2( 



ik-r 



Anr 



d d k 
(2n) d ~ 

<5Cp L (r,T)„ 



1 



k 2 + n 00 (k) kVJ 



(61) 
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The leading-order term in ( 160]) gives a g 6 contribution to SC^^r, T)ni, all other contributions 
being at least of order g 7 /(rT) 2 , 

*C PI Xr,T) m =^ a ^ + yy. (62) 

D. The contribution from diagrams of type IV 

The transverse static-gluon exchange between the two temporal-gluon lines (diagram IV 
and the diagrams derived from IV by inserting gluon self energies in each of the gluon lines) 
receives the following contributions. 

1. Contribution from the scale 1/r 

The contribution from the scale 1/r reads at leading order (with £ = 1) 

,6 AT 2 



SC PL {r,T) wl/t - —-^-C^ J —J —J — e 

(2k! + p) ■ (2k 2 + p) 



k 2 k 2 (k 1 + P ) 2 (k 2 + p) 2 p 2 

Gluon self-energy insertions are suppressed by g 2 . 



(63) 



2. Contribution from the scale T 

The contribution from the scale T vanishes, because scaleless, if no self-energy in- 
sertions are considered. Hence, the leading contribution from the scale T is of order 
g e /T x g 2 T ~ g\ 

3. Contribution from the scale mo 

The contribution from the scale m^ reads 

5C PL (r,T) IVmD = 

# 6 iV 2 -1 , f d d ki f d d k 2 f d d p r „,..., , , x2xl 



4T 2iV 2 

(2k! + p) ■ (2k 2 + p 
(k 2 + m^)(k 2 + m^)((ki + p) 2 + m 2 D )((k 2 + p) 2 + m 2 D )p 



x ^i -r y> ■ k^2 -r V) \l + O(q)} (64) 



This corresponds to a contribution of order g 6 (m£,/T) ~ g 7 , which is beyond our 
accuracy. 
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The leading contribution to <5Cpl(V, T)rv comes, therefore, from £CVl(V, T) IV i/ r , which 
can be computed in dimensional regularization with the help of Eq. ( 156]) . Our final result 
reads 

M*l(t, T) IV = ^^ (l - ^J + O (/) • (65) 



The same result follows from 



20] by expanding in rrriD *C 1. 
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FIG. 6: Diagram V is the first contribution to (Tr^4Q(0)TrAg(r)) c , whereas diagram VI is the first 
contribution to (Tr^(0)TrAg(r)) c . 



E. The contribution from diagrams of type V 

Diagrams contributing to the correlators of six Aq fields in Eq. (I5~TI) are shown in Fig. [61 
The LO diagram contributing to (Trv4Q(0)Trv4Q(r)) is diagram V, which gives 

5C PL (r,T) v - ^ p^. (66) 

If we consider diagram V with gluon self-energy insertions in one of the temporal lines, 
in analogy to ( |59l) . then this starts contributing at order g 7 /(rT) 2 , which is beyond our 
accuracy. 

F. The contribution from diagrams of type VI 

Diagrams contributing to (Tr^^TrA^r)) are like diagram VI in Fig. Eland diagrams 
derived from VI by inserting gluon self energies and other radiative corrections. Colour 
factors aside, their leading contribution may be estimated by simply multiplying the con- 
tribution of the diagrams of Fig. [5] to the Polyakov-loop correlator with the contribution of 
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the diagrams of Fig. [3] to the Polyakov loop. Hence, diagrams of type VI contribute at LO 
to order g A /(rT) 2 x g 2 mo/T ~ g 7 /(rT) 2 , which is beyond our accuracy. 

G. The Polyakov-loop correlator up to order g 6 

Summing up all contributions, we then have 



N 2 y ( rT )2 rT T 

a 3 . N 2 -2 1 a 3 {Z\„ 10 



f (^~6AT- + ^ ^ff [y Ca ~ Y n/ + 27£/3 ° 



+ - 



r y 



& 



C A l-21n^ + 2- — ) H-2n/h,2 

+'«:'^-^7raJC x UC)(v"(//-).-^TT). «w) 



771?, 7T 

where we have made explicit the scale dependence of a s in the leading term. Note that the 
r, T and ma independent term proportional to —27ra 3 CA/9 comes from Eq. (1571) . so it is 
actually a contribution from the scale 1/r that accounts for the matter part of the gluon 
self energy. The term proportional to a 3 /(rT) 3 comes from diagram V, Eq. fl66l) . and from 
the singular part of the gluon self energy in the static gauge, Eq. (J5j 

H. Comparison with the result of Nadkarni 



We compare here with Nadkarni's (N) computation of the Polyakov-loop correlator 20]. 
The regime of validity of Nadkarni's computation is T ^> 1/r ~ mc, while ours is 1/r ^> 
T ^> mo- Therefore, we may only compare results obtained here that do not involve the 
hierarchy rT <C 1, with Nadkarni's results that do not involve the hierarchy rT 3> 1, 
expanded for rrric <C 1. 

In [2o|, the tree-level expression of g 4 (TiA 2 {0)TiA 2 (r)) c /(AT 4 ) reads {N 2 - l)/(8N 2 )a 2 
exp(— 2rm D )/(rT) 2 , which expanded for rm D <C 1 gives SCp^r, T)i, the LO of SCpi,(r, 
T)\im D (to be read from Eq. fl60l) ) and <5Cpl(V, T)ni- Also, the tree-level expression of 
g 6 (Tr74Q(0)TrAQ(r)) c /(36T 6 ) in [20| agrees with SCpi,(r, T) v once expanded for rm D <C 1. 

Diagram IV in Fig. also contributes to Nadkarni's calculation. The diagram does 
not involve gluon self-energy insertions and therefore its calculation does not rely on the 
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hierarchy between 1/r and T. As already remarked, SCpi,(r, T)p/ agrees with Nadkarni's 
result once expanded for rrriD <C l. 8 

Let's now consider the NLO contribution to 5Cpi,(r, T) IImo - This contribution is given 
by the static part of Eq. flBTJl) : 

5C PL (r,T) IINmD --^^-— n J (27r)d(k2 + m|))2 - (68) 

The integral is divergent. In our case, i.e. assuming 1/r ^> T ^> mo, the divergence cancels 
against SCp^r, T)u T , eventually leading to a finite result in 5Cp^(r, T)u T+niD . The In mo/T 
term in Eq. (IBTJj) signals precisely that a divergence at the scale mo has canceled against a 
divergence at the scale T. In Nadkarni's case, i.e. assuming T ^> 1/r ^> mo, we get, along 
with £Cpl(V, T)iiNm D , a contribution from the scale 1/r, which is 

^ PL (r,T) IIN1/r --^^-— /x y^-pe - . (69) 

This is like Eq. (|55|) . but involves only the static part of the self energy (|25|) . since non- 
static modes have been already integrated out at the larger scale T. According to Eq. (|25p . 
we have nj) (|k| ^> mo) ~ T|k| 1_2e . The Fourier transform of l/|k| 3+2e originates a 1/e 
pole. It is this divergence that in Nadkarni's hierarchy cancels against the divergence in 
SCpi J (r,T)n NmD leading to the finite result 



5Cp L (r } T) UNniD +8Cpz J (r,T)nNi/T 



2 1 „,3 



N 2 -la 



3 

ln(2m D r) + 7e - - + 0(rm D ) 



, (70) 



2N rT 

which agrees with the result in [20j. 9 In this case, the lnm^r term signals that a divergence 
at the scale mo has canceled against a divergence at the scale r. 

V. THE POLYAKOV-LOOP CORRELATOR IN AN EFT LANGUAGE 

The calculation of the Polyakov loop discussed in the previous section can be conveniently 
rephrased in an effective field theory (EFT) language that exploits at the Lagrangian level 
the hierarchy of energy scales in Eq. ( l52l) . The EFT framework has the advantage to allow 
more easily for systematic improvements of the calculation and to make more transparent 
its physical meaning. 



8 In Nadkarni's paper this contribution is called ///. 

9 In Nadkarni's paper this contribution is called //. 



23 



Our starting point is QCD with a static quark and a static antiquark, denoted in the 
following as static QCD. Its action in Euclidean space-time reads 

1 IT / n f \ 

Sqcd = J drjd 3 x I ^D ij + X ^D oX + \f« u F* u + ^ q t p qi j , (71) 

where Dq = 8q — igAo, if) is the Pauli spinor field that annihilates a static quark, \ is the 
Pauli spinor field that creates a static antiquark, and q±, ..., q nf are the light quark fields, 
which are assumed to be massless in this study. 

The Polyakov-loop correlator may be expressed in static QCD as 

(Tr4(0)frMr)> = ^<x](0, 1/T)^(r, 1/T)$(r, 0)^(0, 0)), (72) 

where M = [<5 3 (0)] 2 and we have written explicitly the colour indices. The thermal average 



on the right-hand side reduces to the Polyakov-loop correlator on the left-hand side a: 

n ft 

integrating out the fields ip and x [2|- On general grounds, one also expects that |38| 



ter 



(Tr4(0)frL F (r)) = ± £ e"^, (73) 

n 

where E n are the eigenvalues of the Hamiltonian associated to the static QCD Lagrangian. 

A. pNRQCD 

Potential non-relativistic QCD (pNRQCD) is the EFT that follows from QCD by in- 
tegrating out from the static quark- antiquark sector gluons of energy or momentum that 
scale like the inverse of the distance r between the quark and the antiquark. Since 1/r is 
the largest scale, the matching of the pNRQCD Lagrangian may be done by setting to zero 
all other scales and, in particular, the thermal ones; as a consequence, the Lagrangian is 



identical to the one derived at zero temperature 34], l40l-l43| . In Euclidean space-time, the 
action reads 

SpNrqcd = f dr J d 3 x j d 3 r Tr J S\d + V S )S + O\D + V o )0 

-iV A (S f r ■ #EO + O f r • #ES) - l -V B (O f r • #EO + O f Or ■ gE) 

+ % Vc (rW&gEPO - rVO+OD^^') + ^ pNRQCD 
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1/T 



dr 



n f 



d 3 x I --F« v f; v + ^J gtfto 



(74) 



where the trace is over the colour indices, S = 1n x nS/vN is a quark-antiquark field in 
a colour-singlet configuration, O = \/2T a O a is a quark-antiquark field in a colour-octet 
configuration, D O = <9 — ig[A ,O], D = V — igA and E l = F i0 is the chromoelectric 
field. The fields S and O a depend on the continuous parameter r that labels the distance 
between the quark and the antiquark, the centre-of-mass coordinate x and the Euclidean 
time r; the gluon fields have been multipole expanded and, therefore, depend on x and 
r only. The quantities V s , V , Va, Vb and Vc are the matching coefficients of the EFT. 
These are non-analytic functions of r. Since Va(j') = 1 + 0(a 2 ) J44|, Vb[t) = 1 + 0(a 2 ) 
and Vcir) = 1 + 0(a s ) it will suffice to our purposes to put Va{t) = Vsir) = Vc(r) = 1 
from now on. VI and V are the singlet and octet potentials in pNRQCD: V s is known up 



to three loops 



45] and V is known up to two loops |46| . For the purpose of obtaining the 



Polyakov-loop correlator at NNLO accuracy it is sufficient to know V s and V at one-loop 
accuracy and their difference at two-loop accuracy: 



V.(r) 
V (r) 



-a 



« s (l/r) 



1 



31 
~9 



Ca ~ Y n f + 2 7^o) ^ + 0(a 



1 a s (l/r) 

2N r 



l+[*CA-^f + *r E fh)% + 0{a 



(N 2 - l)K(r) + V a (r) 



N(N 2 -l)a; 



TV' 



- 3 [1 + 0(a B )] 



(75) 
(76) 
(77) 



/■ v 4 

Finally, <5£ p nrqcd includes all operators that are of order r 3 or smaller. At tree-level, they 
may be read from the multipole expansion of the quark and antiquark coupling to the tem- 
poral gluon in the static QCD Lagrangian ( I7ip . hence they just involve covariant derivatives 
acting on a chromoelectric field: the leading-order operator being —ir l r^r k TT{O^D l D : >gE k S 
+S jf D l D : 'gE k O}/24: 42j. As we will argue in the next section, these terms are of order g A , 
however, their contribution eventually cancels in the Polyakov-loop correlator up to order 
g 6 (rT)°. For this reason, we do not need to specify them further here. 
Matching the connected Polyakov-loop correlator to pNRQCD gives 



C PL (r,T) 



1 

iV2 



_ (S(r,0,l/r)fft( r ,o,Q)) (Q°(r, , l/T)O a \v, 0, 0)) 
■^s — r £ 



M 



N 



+0 (« s 3 (rT) 4 ) 



- (L F ) 2 . 



(78) 
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The right-hand side is the pNRQCD part of the matching. It contains the singlet and octet 
correlators, (5(r, 0, l/T)S^(r, 0, 0)) and (O a (r, 0, l/T)O a ^(r, 0, 0)), not surprisingly because 
in the r — > limit the tensor fields x](0, l/T)ipi(r, 1/T) and ip\ [(r, 0)xj(0, 0), appearing in 
the right-hand side of Eq. (1721) . decompose into the direct sum of a colour- singlet and a 
colour-octet component. The colour-singlet and colour-octet correlators may be read from 
the Lagrangian (1731) : 

(O-(r,Q,l/r)O-t(r,0,0)) = e . v . (rWT [(Jy2 _ 1} (La) + ^ j (g0) 

where £ s and <5 stand for loop corrections to the singlet and octet correlators respectively. 
The factor (La) comes from the covariant derivative D acting on the octet field in f J74|) . 10 
Note that at finite temperature, for T > g 2 /r, the octet correlator is not suppressed with 
respect to the singlet one, while in the opposite limit, T <C g 2 /r, the Polyakov-loop corre- 
lator is dominated by the singlet contribution. Higher-dimensional operators have not been 
displayed, because they are negligible with respect to our present accuracy, which is of order 
g 6 (rT)°. The reason is that higher- dimensional operators involve the coupling with at least 
two field-strength tensors, hence the corresponding matrix elements are at least of order 
(rT) A ; moreover, as can be seen by adding two external gluons to diagram I of Fig. [51 the 
matrix element of an operator coupled with two external gluons is at least of order g 6 . The 
normalization factors Z s and Z have to be determined from the matching condition (1751) . 
While V s and V are the same at zero and finite temperature, the normalization factors are 
not for they depend on the boundary conditions. 

In order to determine the normalization factors Z s and Z , let us consider in Eq. ( 178]) only 
contributions coming from the scale 1/r. In dimensional regularization, all loop corrections 
vanish in the pNRQCD part of the matching and the Polyakov loops (Lp) and (La) reduce 
to one; therefore, the matching condition reads 

C PL (r,T) 1/r = (Tr4(0)frL F (r)) 1/r -l = ± [Z s e- y ^ T + Z (N 2 - l)e^W/ T ]-l. (81) 



10 The adjoint Polyakov loop (La) factorizes the contribution coming from the gluons in the thermal bath 
that bind with the colour-octet quark- antiquark states to form part of the spectrum appearing in the 
right-hand side of Eq. (|T3[) . In pNRQCD at zero temperature, a similar expression factorizes the non- 



perturbative gluonic contribution to the gluelumps masses 34 [ 

2(3 



We may now proceed in different ways. A way consists in matching with the spectral 
decomposition f!73|) . By noting that at the scale 1/r the spectrum is just given by a singlet 
state of energy V s (r) and iV 2 — 1 degenerate octet states of energy V (r), the matching 
condition implies that Z s — Z — 1. Another way consists in taking advantage of the 
Polyakov-loop correlator calculation done in Sec. IIVI and matching to it. Cpl(V, T)i/ r is the 
sum of Eq. (!53j) . Eq. (l57|l without the contribution from the matter part of the gluon self 
energy, Eq. fl65l) and Eq. fl66l) ; it reads 

_ , N 2 -l (a s (l/r) 2 al N 2 -2 1 « b 3 /31_ 10 

4 c *H)+°{M- (82) 

A direct inspection shows that this expression satisfies 

C PL (r, T) 1/r = ± [e~ v ^ T + (N 2 - l)e~ v °^ T ] - 1, (83) 

up to order a%, for V s (r) and V {r) given by Eqs. (f75l) - (jTTj) . n We note that Eqs. (1821) and 
( |83l) are equivalent for rT ^> (yf 2 , however, in Eq. ( !8~3l) . we resum some contributions that 
would become large for rT < g 2 . Equation (1831) is therefore valid also in that regime. Finally, 
we observe that the combination of the two procedures provides a non-trivial verification of 
Eq. (I77j) . i.e. of the two-loop difference between the octet and the singlet potentials, known, 
so far, only from the direct calculation of the two-loop octet potential in a covariant gauge, 
done in Ref. [46|] . 

Loop corrections to the singlet and octet correlators in Eqs. ( |79l) and (IHDl get contribu- 
tions from the scales T, mjj and lower ones. We now proceed to evaluate these corrections, 
separating the contributions of the temperature from the ones of the Debye mass. 

B. The temperature scale 

In the hierarchy ( f52|) . the next scale after the inverse distance is the temperature. Our 
aim is thus to compute the temperature contributions to loop corrections in pNRQCD. These 
loop corrections are the terms S s and 5 that were introduced in Eqs. (!79|) and (IHUj) . We call 



11 More precisely, the matching to (|52|) fixes Z s = Z = 1 up to order o^ and Z s + (^V 2 — 1)Z = N 2 up to 
order a^. 
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5 St x and 8 ,t the parts of 5 S and 5 respectively that encode the contributions coming from 
the scale T; they may be obtained by expanding 5 S and 5 in mo, V s , V and in any lower 
energy scale. Similarly, 5(L R ) T is the part of (L R ) that encodes the contributions coming 
from the scale T. Different terms contribute to S Sj t, <5 ,t and 5(Lr)t] we examine them in 
the following. 




FIG. 7: The pNRQCD Feynman diagram giving the leading-order correction to 5 S . The single 
continuous line stands for a singlet propagator, the double line for an octet propagator, the circle 
with a cross for the chromoelectric dipole vertex proportional to Va in the Lagrangian (|74p and 
the curly line connecting the two circles with a cross for a chromoelectric correlator. 

1. The singlet r 2 contributions 

We start considering the onedoop, order r 2 in the multipole expansion, correction to 
the singlet correlator induced by the diagram shown in Fig. [7J it reads 




x 



e _ i(T _ r>n {E ia UabE j b) ^ k y (g4) 



In the sum integral, we may distinguish between contributions coming from the non- 
zero modes and from the zero modes. 

For the contribution coming from the non-zero modes, only the leading-order chromo- 
electric correlator in momentum space (E ia U a bE^ b )(u n ,k) (U a b stands for a Wilson 
straight line in the adjoint representation connecting E ia with E^ ; at leading order 
U a b = o~ab) is relevant at our accuracy: 



(E ia U ab Ei b )(uj n ,k) = (N 2 -l) 



+ (Sij - k l k 3 )- 



2 
" n 



(85) 



k 2 v lJ y w 2 +k 2 _ 

Loop corrections to the chromoelectric correlator contribute to the Polyakov-loop cor- 
relator at order g 6 (rT) or smaller. Because of the hierarchy (|52|) . we can expand 
the right-hand side of flHll) in V — V s . The longitudinal part of the chromoelectric 
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correlator, i.e. the first term in square brackets, vanishes in dimensional regulariza- 
tion, whereas the transverse part is sensitive to the scale T through the Matsubara 
frequencies. After performing the sum integral over the non-zero modes, we obtain 

*S r2)NS = -9 2 c/^(V -V s )+g 2 C F ^(V -V s r + o(g°(rT),^ 

= - 2 KNC F a 2 s rT + ^N 2 C F a 3 s + O (g\rT), ^-] . (86) 



9 s 36 s V rT , 

The contribution coming from the zero modes reads 

'#"' = ( !9 ^) ^/(gl(^"C^'>(0.k)ll^r + O( 9 6 (rT)). (87) 

Here, the first no n- vanishing contribution in dimensional regular ization comes from 
the one- loop correction to the chromoelectric correlator. The integral with (E ia 
U a bEi b )(0, k) at one loop has been calculated in [9J]. Using that result we obtain 12 

Sff )S = \aZ)C F ^{r mD y - 2 -amC F al{rTf + O (g\rT)) . (88) 

2. Higher multipole terms 

Our aim is to calculate in the EFT the Polyakov-loop correlator at order g 6 , neglecting 
terms of order g 6 (rT) or smaller. Contributions coming from the <5£ p nrqcd P ar t of 
the pNRQCD Lagrangian, which includes terms of order r 3 or smaller coming from 
the multipole expansion, share, at leading order, the same colour structure and the 
same order in a s as Eqs. flHBT) and (188]) but are suppressed by powers of rT. We may 
write these contributions as 

oo oo 

^ pN rqcd = ^ 2 )NS^ c NS (rT) 2n +2 + ^)S^ c S (rT) 2n + 2 + ^6 (rT) 3) 5 

n=0 n=0 



12 The chromoelectric correlator is gauge invariant. In static gauge, thermal corrections arise from the 
non-static part of the spatial gluon propagator. Hence, at one loop, only gluon self-energy diagrams may 
provide thermal corrections; we have 

k^ +ng l ( f(k) mat 

where n^Q S (k) ma t is the matter part of the gluon self-energy's temporal component calculated in static 
gauge, which can be read from Eq. (flT]) . Finally, we recall that n^ s (k) mat is the same in static gauge 
and in Coulomb gauge. 
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(89) 

where the unknown coefficients c^ s and c^ are, as we will see, irrelevant for the purpose 
of calculating the Polyakov-loop correlator at order g 6 (rT)°. 

3. The octet contributions 

As in the singlet case, one loop-corrections to the octet correlator may be divided into 



order r 2 non-zero mode contributions (<5 



0(r 2 ) NSx 
o,T J 



order r 2 zero-mode contributions 



(K 



0(r 2 )S 



i, and higher multipole terms (<5 oT pNRQCD ). It turns out that 



,0(r 2 ) NS 
>o,T 



s 



C(r 2 )NS, 



\Vs*r^Vo1 



(90) 



and, up to order g 6 (rT)°, 



S 



0(r 2 )S 
o,T 



S 



0(r 2 )S 
s,T ' 



r <5£pNRQCD 

1 o,T 



-5 



<5-CpNRQCD 



(91) 
(92) 



These equalities are proved in appendix | 



4. S(L R ) T 

Finally, we need to calculate the contributions to the Polyakov loop coming from the 
scale T. The order g A contribution may be read from Eq. ( 139]) . Since we do not know 
the order C^g 6 contribution, we write o~(Lr)t as 



S{L R ) 



CrOl 2 s 



T 



C A In 

2e 



4T 2 
1* 



1 — 7b + ln(47r) J — nf In 2 + a a 



+ O (at) , 
(93) 



where the explicit value of the coefficient a does not matter. Instead, what matters 
here is that this coefficient is common to all colour representations. The first correction 
from the scale T not of the type Cr a™ appears at order a^ and comes from diagram b) 
in Fig. H]with two self-energy insertions, one in each temporal gluon. Note that Eq. ( 14H1) 
provides the first correction not of the type C^a™ coming from the scale mo- 

In summary, we obtain the contribution of the scale T to the singlet and octet correlators: 

e -Vs(r)/T 6sT = 



-Vs(r)/T < 



- -nNC F a 2 s rT 
9 



i + £ c « NS (^) 



2n+2 



n=0 



TV 
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N 2 C F al 
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( h{Z)C F -{rm D f - \amC F al{rTf 
L ix 3 



l + £c£(rT) 



2n+2 



n=0 



+0(g b (rn^ 



(94) 



e -Vo(r)/T^ N 2_ 1)S{LA)T + So ^ 
{N *- 1)e -V (r)/ T iCA a 2 



1 4T 2 

C A ( g- - In —j- + 1 - 7b + ln(47r) ) - n f In 2 + a a s 



+ -7ra„rT 
9 



i + E c « s (- T ) 



2n+2 



n=0 



+ f 2 Na! 



l((^-(rm D y - k(3)« s 2 (rT) 2 
4 TV 7r 3 



l + ^4(rT) 



2n+2 



n=0 



+0U 6 (rT), 



r y 



(95) 



Inserting Eqs. ()93j)-(l95l) into Eq. ( 1781) and expanding, we obtain that the connected 
Polyakov-loop correlator is given by 



C PL (r,T) = C PL (r,T) 1 



It 



1 ai 



*_ 3 W 

'18 F " s+ 8iV2 rT 



C A 



1 u 2 
2 In Ar 

e 4T 2 

+2n/ln2 



2 + 27^-2 ln(47r; 



^K(rT), 



w 



+ loop corrections at the scale mo or lower 



(96) 



where Cp^r, T)y r may be read from Eq. ( 182|) . We observe that, in the connected Polyakov- 
loop correlator, terms proportional to the unknown coefficients c„ , c„ and a have canceled. 
The thermal corrections in ( 1961) agree with those calculated in Sec. IIV( in particular, they 
correspond to the sum of the gluon self-energy matter-part contribution in Eq. ( |5"7|) with 
Eq. (!58|) . The result in Eq. (!96|) has an infrared divergence that originates at the scale T. 
This divergence shall cancel against an opposite ultraviolet one at the scale m^, which will 
be the subject of the next section. 



C. The Debye mass scale 

Here we compute the contributions to the singlet correlator, the octet correlator and the 
Polyakov loop coming from loop momenta sensitive to the Debye mass scale. We call these 
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contributions S s ^ mD , 0~ OtTnD and 5{L R ) mD respectively. They may be computed by evaluating 
the loop integrals in 8 S , 5 and 5(L R ) over momenta of the order mo and expanding with 
respect to any other scale. The Debye mass scale is the lowest scale we need to consider 
here; contributions coming from scales lower than mu are beyond our accuracy. Different 
terms contribute to 8 s ^ mD , S 0jmD and 5{L R ) mD ; we examine them in the following. 

1. The singlet and octet contributions 

The leading-order contribution to <5 s , mD comes from the self-energy diagram shown in 
Fig. [7] when evaluated over loop momenta of order mo- The contribution reads 

xe -i(r-r')u> n (E ia U ab Ei b )(u n ,k)\\khm D - (97) 

The chromoelectric correlator evaluated over the region |k| ~ mo gives rise to scaleless 
momentum integrals unless for the temporal part of the zero mode, n = 0, which is at 
leading order (E ia U ab E^ b )(0, k)|| khmD = (iV 2 - 1) k { k j /(k 2 + m 2 D ). We obtain 

(98) 



J s,m D 



The leading-order contribution to 0~ O)TnD comes from the octet self-energy diagrams 
shown in Fig. |9j when evaluated over the region |k| ~ mo- Also in this case, the 
only non-vanishing contribution comes from the zero mode of the temporal gluon 
propagator, which is l/(k 2 + m 2 D ) (see Eq. (122]) ). For the same argument developed 
in appendix [Gl we find that 

TY1 - 

Higher multipole terms are of order a s r 2 [rmo) 2 ~ g 7 (rT) 4 or smaller and, there- 
fore, beyond our accuracy. 

2- o~(Lft) mD 

We need to calculate the contribution to the Polyakov loop coming from the scale mo ■ 
It may be read from Eqs. (1331) . (jSJ) and (jUD- Since we do not know the order C^g 5 
and C^g 6 contributions, we write (LR) mD as 
C R a s m D 



8{L R ) mD 



T 
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C R a\ 
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Cxl-^-i"£r-5 + 7--M^ 



big + b 2 g 2 



+ ( 3 c|-^i)|(^) ! + o( 9 T, 



(100) 



where the explicit values of the coefficients b% and 62 do not matter. Instead, what 
matters here is that these coefficients are common to all colour representations. 

In summary, we obtain the contribution of the scale mo to the singlet and octet correlators: 

(101) 



e -W T 5 



-V s i 



s,m,E> 



e-^)lA-C F ^r^ + 0(g\rT))\ 



-V.M/T [{N 2 _ 1} 6{LA)mD + ^^ = (N 2 _ 1)e -,( r )/T|^| + ^ clal ^y 
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Vt J 



C A a 



C A In 

2e 



£^--- + jE-H^))+i>K/ + i>2<f 



"D 



+l^r 2 ^ + o(, 7 ) y 



(102) 



Inserting Eqs. f ll00p - (U02p into Eq. ([55]) and expanding, 13 we obtain that the connected 
Polyakov-loop correlator is given by 

C PL (r,T) = C PL (r,T) 1/r 

- 18 ^ + ^iV 2 "" 8 IT") 



+ 



N 2 -la. 



N 2 rT I 4 T 
,7 



T 



T 2 



^(--ln— + -)-n / ln2 



m 2 ' 2 



+0 ^(rT), 



£ 



(103) 



( rT )2 / 

where CpL(r,T)i/ r may be read from Eq. (|82|) . We observe that, in the Polyakov-loop cor- 
relator, terms proportional to the unknown coefficients b\ and 62, as well as the divergences, 
have canceled. The origin of the thermal corrections to the Polyakov-loop correlator in the 



13 In terms of 8 8>T , 5 s ,m D , 5 q ,t, $o,m D , 5{L F ) Tl S(L F ) mD , S(L A ) T and S(L A ) mD , Cp L (r,T) reads 



+e -v (r)/T ^ N 2 _i)(i + s(L A ) T + 5{L A ) mD ) + S , T + 5 . mD ] } 



(l+5(L F ) T + 6(L F ) mD ) 2 
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situation 1/r 3> T ^> m£> ^ g 2 / 7 ' is clear. The term — CV7rag/18 arises from the dipole 
interaction contributions and from their interference with the zero-temperature potentials. 
The other thermal corrections arise from the interference of the adjoint Polyakov loop with 
the zero-temperature potentials. 

The result coincides with Eq. (16"Tj) . obtained in Sec. [TV] after a direct calculation. The 
differences in the way the two results were achieved illustrate well the typical differences 
between a direct computation and a computation in an EFT framework. In the EFT frame- 
work, some more conceptual work was necessary in order to identify the relevant contri- 
butions. Once this was done, we could take advantage of previously done calculations (in 
particular for V s (r) and V (r)) and reduce the calculation to essentially one diagram, shown 
in Fig. evaluated in different momentum regions. In the EFT framework, we could also 
gain some new insight by reconstructing the spectral decomposition of the Polyakov-loop 
correlator and by providing two new quantities: the colour-singlet and the colour-octet 
quark-antiquark correlators. 

D. Singlet and octet free energies 

Potential NRQCD at finite temperature allows to define a colour-singlet correlator, 
(S(r, 0, 1/T)5 f (r, 0, 0)), and a colour-octet correlator, (O a (r, 0, 1/T)0 at (r, 0, 0)), which are 
both gauge-invariant quantities. We may associate to them a colour- singlet free energy, 
f s {r,T,m D ), and a colour- octet free energy, f {r,T,m D ), such that 



{Sir, 0,l/T)St( r ,0,0)) = e- v °W T (l + 5 s>T + 5 SimD ) 



= e ~fs(r,T,m D )/T 



(104) 
(O a (r, 0, l/T)O at (r, 0, 0)) = e "^)/ T [(N 2 - 1) (1 + 6(L A ) T + 8{L A ) mD ) + 5 , T + 5 , mD ] 

= (N 2 - l) e -Mr,T,m D )/T_ (1Q5) 

Using the results of the previous sections, we have that 



f s (r,T,m D ) = V a (r) 



+-nNC F a 2 rT 2 
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as „2^3 | T /n / „G/^rr\ 9 
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(107) 



We note that / s (r, T, mc) and / (r, T, mo) are both finite and gauge invariant. They also 

do not depend on some special choice of Wilson lines connecting the initial and final quark 

and antiquark states. 
f] 
In [9[, the colour- singlet quark-antiquark potential was calculated in real-time formalism 

in the same thermodynamical situation considered here and specified by Eq. (I52p . The result 



may be found in Eq. (92) of Ref. 9j. Comparing terms of the same order, the real part of the 

1 71 

real-time potential differs from f s (r, T, mo) by -7rA^Cj?Q; 2 rT 2 N 2 Cfo^T. The origin of 

the difference may be traced back to terms in Eq. ( 1M|) that would vanish for large real times. 
Indeed, performing the calculation of (5*(r, 0,t)S'(t, 0, 0)) for an imaginary time r < 1/T, 
along the lines of Sees. IV Bl and IV CI and then continuing analytically r to large real times, 
one gets back exactly both the real and the imaginary parts of the real-time colour- singlet 
potential derived in J9J. The difference between the singlet free energy and the real part of 
the real-time colour- singlet potential appears to be a relevant finding to be considered when 
using free-energy lattice data for the quarkonium in media phenomenology. 



E. Comparison with the literature 



An EFT approach for the calculation of the correlator of Polyakov loops was developed 
in 47J for the situation mr> > 1/r and in 20j for T ^> 1/r. In neither of the two cases, 
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the scale 1/r was integrated out: the Polyakov-loop correlator was described in terms of 
dimensionally reduced effective field theories of QCD, while the complexity of the bound- 



state dynamics remained implicit in the correlator. The description developed in [20|, |47| is 
valid for largely separated Polyakov loops. Under that condition, the correlator turns out 
to be screened either by the Debye mass, for riri£, ~ 1, or by the mass of the lowest-lying 
glueball, for rmc 3> 1. 



In [38| , the spectral decomposition of the Polyakov-loop correlator was analyzed. It was 
concluded that the quark-antiquark component of an allowed intermediate state, i.e. a 
field if describing a quark located in xx and an antiquark located in x 2 , should transform 
as </?(xi,x 2 ) — > g(x 1 )y9(x 1 ,x 2 )(^(x2) under a gauge transformation g. Equation (175|) is in 
accordance with that result for, in pNRQCD, both the singlet field S and the octet field 



O transform in that way 4Gj. We remark, however, a difference in language: in our work, 



singlet and octet refer to the gauge transformation properties of the quark-antiquark fields, 



while, in [38|, they refer to the gauge transformation properties of the physical states. 

In 31] , a weak-coupling calculation of the untraced Polyakov-loop correlator in Coulomb 
gauge and of the cyclic Wilson loop was performed up to order g 4 . Each of these objects 
contributes to the correlator of two Polyakov loops through a Fierz transformation that also 
generates some octet counterparts. It is expected that large cancellations occur between 
those correlators and their octet counterparts in order to reproduce the Polyakov-loop cor- 
relator given in Eq. (!67l) . Such large cancellations should occur at the level of the scales 
1/r, T and nip as we have already experienced in this work. Note that in the case of the 
untraced Polyakov-loop correlator, the octet contribution shall also restore gauge invariance. 

VI. SUMMARY AND OUTLOOK 

In the weak-coupling regime, we have calculated the Polyakov loop up to order g A and 
the correlator of two Polyakov loops up to order g 6 (rT)°, assuming the hierarchy of scales 

i „, g 2 

- > T>m D > — . 



The Polyakov-loop calculation differs from the result of Gava and Jengo [2l| by a finite 
contribution at order g A . We have analyzed in detail the origin of the difference and shown in 
an appendix that our result may be reproduced also performing the calculation in Feynman 



gauge. Our calculation agrees with the recent finding of Ref. [31]. 
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The calculation of the Polyakov-loop correlator is new in the considered regime, although 
some partial results may be deduced from a previous work of Nadkarni, who studied distances 



r ~ I/mo [201 ] . We have performed the calculation in two different approaches: by a direct 
computation in static gauge and by calculating the Polyakov-loop correlator in a suitable 
EFT that exploits the hierarchy of scales in the problem. In this second approach, we 
have used pNRQCD at finite temperature and subsequently integrated out lower momentum 
regions. The advantages of this second approach are that the calculations do not rely on any 
specific choice of gauge and the systematics is clearer. Moreover, it makes explicit the quark- 
antiquark colour-singlet and colour-octet contributions to the Polyakov-loop correlator. In 
particular, we have shown that at leading order in the multipole expansion the Polyakov-loop 
correlator can be written as the colour average of a colour- singlet correlator, which defines 
a gauge-invariant colour-singlet free energy, and a colour-octet correlator, which defines a 
gauge-invariant colour-octet free energy. This is in line with some early intuitive arguments 



given 



m y,y, 



201 ] . In general, however, such a decomposition does not hold. 

In the weak-coupling regime, the degrees of freedom of pNRQCD are quark-antiquark 
colour- singlet fields, quark-antiquark colour-octet fields, gluons and light quarks. The ob- 
tained result for the Polyakov-loop correlator is consistent with its spectral decomposition. 
In the strong-coupling regime, the degrees of freedom are expected to change when the typi- 
cal energy of the bound state is smaller than the confinement scale Aqcd- in that situation, 
the bound state would become sensitive to confinement and give rise to a new spectrum 
of gluonic excitations (hybrids, glueballs). In the present work, we have not discussed this 
situation, which surely deserves investigation. 

Possible further extensions of this work also include the study of the Polyakov-loop cor- 
relator in different scale hierarchies, in particular at temperatures of the same order as or 
ligher than 1/r, where the present analysis should smoothly go over the ones performed in 
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471 ] . As mentioned above, also analyses that involve the strong-coupling scale should be 
addressed. 

Finally, the present study should be completed by the study of correlators different from 
the Polyakov-loop one. Among these, the most studied in lattice gauge theories are the 
untraced Polyakov-loop correlator and the cyclic Wilson loop. Also the octet Wilson loop 
should be included for its role in the Polyakov-loop correlator. Since some partial perturba- 
tive results are already available for some of these correlators, it would be interesting to see 
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how they can be reproduced in the EFT framework introduced here and how they combine 
to give back the Polyakov-loop correlator. 
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Appendix A: Feynman rules in the static gauge 

In the following, we list the Feynman rules in Euclidean space-time under the gauge 
condition do A = 0. The temporal propagator reads (dropping colour indices) 

A,oK,k)= = |f, (Al) 

where, as usual, u n = 2imT and the Kronecker delta fixes n = 0, making this propagator 
purely static. The spatial propagator can be divided into a non-static {n ^ 0) and a static 
(n = 0) part. The former reads 

1 / kk\ 

D tJ (u n ^ 0, k) = k£&QQQQQS£L> = 2 M2 6 v + ^T H 1 ~ <W> ( A2 ) 

and thus mixes longitudinal and transverse components. The static part has a residual gauge 
dependence on the parameter £; it reads 

A>n = 0, k) = <-^~-v~s~^ = -L U 3 _ (i _ £)^\ S n0 . (A3) 
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Finally the ghost propagator reads 

-Dghost(w n , k) 



<5n0 
k 2 " 



(A4) 



and is thus purely static 14 . The interaction vertices (gluon-gluon and gluon-ghost) are the 
usual ones. 



Appendix B: The gluon self energy in the static gauge 

We proceed to the computation of the Matsubara sums in Eq. ((Tj) in order to obtain Eqs. 
ffTUj) . ffTTp . (JT2D and ([13]). We recall the two basic bosonic Matsubara sums 29] 



r£ 



1 



l + 2n B (|p|) 



P 2 +CJ2 



2|p| 



(Bl) 



+oo 



^E 



(p2 + W*)(q» + W*) 2 | P ||q| 



l + n B (|p|) + n B (|q|) , n B (|q|) - n B (|p|) 



IpI + |q| 



|p| - q 



(B2) 



where n B is the Bose-Einstein distribution. Since the sums include also the zero mode, in 
evaluating the master sum integrals defined in Eqs. (J7J) and (jHJ) we will have to subtract it. 
Furthermore, we identify the temperature-independent part (the unity) in the numerators 
on the r.h.s of Eqs. fIBll) and flB2j) as the vacuum part and the part proportional to the 
thermal distributions as the matter part. 
For I , we have 

J P V 



,2 » 



2<- 



d d p fl + 2n B (\p\) T_\ T^ 

,2 J - 12 ; 



(B3) 



(2n) d \ 2|p| p- 

the subtracted zero mode along with the vacuum part vanish in dimensional regularization. 
For ii, we have (q = k — p) 

d d p 



h 



/' 



2c 



p| ^l + ra B (|p|)+7i B (|q|) | n B (|q|)-n B (|p lA T 



IPI + q 



IPI - q 



,2c 



f 



P^ 



IpI^b(IpI) / -2|q 



(27r) d [2|q 
d d p 

W? L2|p||q|(|p| + |q|) ' 2|q| Vp 2 -q 2 

, |q>B(|p|) ( -2|q' 



q^ 



2|p| 



p2 _ q/2 



q 



14 



The non-static ghost can be shown to decouple [22 1 
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where we have operated a shift p— )■ q' = p + k, q— > — p in some terms of the matter part. 
The vacuum part can be brought into a more standard form by noting that 

f +00 dpo 1 _ 1 

;_«, 27r(p 2 +p 2 1 )(q 2 +p 2 ) 2 |p||q|(|p| + |q|)- l J 

This allows to write the three-dimensional integral as a standard Euclidean four-dimensional 
integral, which can be computed with the formulas listed in appendix IB II setting d + 1 = 
4 — 2e. We thus have 



{h) vac — H 



2c 



d d+1 p pVOW - 8^8, 



vO) 



p2g2 



(V - <Vo<Wa* e Lf +1 (k, 1, l)| fe o =0 . (B5) 



(27r) d+1 

The zero-mode integral vanishes in dimensional regularization, whereas the remaining matter 
part is finite and gives 



(/ 



1 J mat 



2tt 2 



d\v\ IpIMIpI) 1 



iPl 
2|k| 



In 



|k| + 2|p| 



|k| -2|p| 



(B6) 



Analogously, we have for I 2 
d d p 



J, _ 1,2.. 26 



The vacuum part is 



I fl + n B (\p\)+n B (\q\) | n B (|q|)-n B (|p lA T 



2|p||q 



ipi + q 



ipi - q 



p2q2 



(^2jvac — H 



2c 



d d+l p k 2 

(27r) d + 1 pV 



k /i e L d+ i(fc,l,l)| fc0=0 , 



the matter part is 



1 / f°° Ikl 

(h)mnt = 2^2 ( / ^IpI ^b(|p|) — In 



|k|+2|p| 



|k| -2|p| 



and the subtracted zero-mode part is 



d: 



2 /zero 



/i 



2r 



d d p Tk 2 



(27r) d p 2 q 2 ' 

which has been kept in dimensional regularization. 
We consider now I^: 

d d p [l + 2n B (|p|) T 



(B7) 



(B8) 



(B9) 
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h = kV 

C^vac = kV 
(-<3,)mat = 



{2ll) d 

<i d p 1 



2|p| 3 
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P J 
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— / d|p||p|n B (|p|) — • 
tH Jo P 



(BIO) 
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(B12) 
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In dimensional regularization the subtracted zero mode vanishes. The matter part is infrared 
divergent. Since this divergence will cancel against terms from I a in the sum fl6]), we present 
the result directly in the three-dimensional limit. 
Ia is given by 

k 4 r k 4 r k 4 



j _ ja jb jc 
1a — 1a — 1a — 1, 



! " ' ! ' J \ pV^ 



p 2 q 2 p 2 



P p 2 qV 
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96T' 
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which is a term peculiar to this gauge; it is singular in the T — > limit and constitutes 
TT NS (k1 • T b is 

ii 00 Wsmg- J4 io 

1 A + n B (|p|)+n B (|q|) , n B (|q|)-n B (|p|)\ T 



4,.2e 



k 4 /i 
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The vacuum part can be brought into a more familiar form by adding and subtracting 
l/(2|p|V) 



\l^)vac 
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Although the matter part of I a is infrared divergent, its infrared divergence cancels against 
the matter part of I3, i.e. Eq. ( 1B12I) . in the sum ((6]). Hence, we may evaluate it directly 
in three dimensions. In contrast, we will keep regularized the subtracted zero modes. As 
discussed in the main text, these subtracted zero modes behave like e|k| 1_2e and are going 
to contribute when evaluating the Fourier transform of |k| 1-2<: /|k| 4 in the Polyakov-loop 
correlator calculation, like in Eq. ( 155|) . Therefore, (/^mat and (iDzero rea d 



i-^Jmat — 


2vr 2 


(, 1 4 ) zero 
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1 f°° Ikl 3 
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Similarly J| reads 
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(B19) 
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1/4 J mat 
1/4 J zero 



4.,2e 



kV 



2tt 2 



d d p 
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(27r) d 2|p| 2 |q| 3 

Ikl 3 

d|p||p|M|p|)7Ji-i3 ln 

2|p| 
d d p Tk 4 



/i ie L d (k, 1,3/2), 
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(B20) 
(B21) 
(B22) 



(27r) d p 2 q 4 ' 

Notice that, as we anticipated, the sum (/3) ma t/2 — (iDniat/4 — (l£)ma,t/4, which is the 
combination appearing in Il^ s (k), is infrared finite. It is also worthwhile noticing that the 
vacuum parts (i|) va c and (I|) V ac are infrared divergent, but that in the sum (/ 3 ) vac /2 — 
(it) vac/4 — (/|)vac/4, these infrared divergences are canceled and replaced by an ultraviolet 
divergence eventually removed by renormalization. The canceling infrared divergence and 
the remaining ultraviolet one come from (i3) vac , which vanishes, like in Eq. (IBllj) . if the 
two are set equal, as usually done in dimensional regularization. 
Putting all pieces together in Eq. (jfJJ) and using 

d d p k 4 



/' 
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(27r) rf p 4 q 2 



Ik^-V 6 ^) 



-3/2 + e r(3/2 + e)r(l/2-e)r(-l/2-e) 
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Ikp-V* 
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/x 2e (47r 



. 3/2+e r(l/2 + e)r(l/2-e) 



Ikl 1 " 2 ^ 



r(l-2e) 
[1 + e(-7B + 1h(16tt)) + C(e 2 



we obtain Eqs. (JTUj) . (ITT]) , (1T2]) and <jT3j> . 



1. One-loop integrals 

We list here the loop integrals Ld, L^ and L^ u , obtained with the Gegenbauer polynomials 
technique J48| : 

d d p 1 



L d (k,r,s) 



L%(h,r,s) 



(2ir) d (p + k) 2r p 2s 
k d-2(r+s) r ( r + s _ d/2) T (d/2 -s)T (d/2 - r) 
(47r) d / 2 r(r)r(s) 



T{d — s — r) 



-W 



(2vr) d (p + k) 2r p 2s 
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(47r) d / 2 r(r)r(s) 
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d [ ,T,S) ~ J (2ny( P + k)^s 
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(47r) d / 2 
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(B25) 
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Appendix C: Expansions 

In this appendix, we list the expansions of the gluon self energy for temperatures much 
greater or smaller than the momentum k. 

We start with T ^> |k|. In the non-static sector, Jo gives its exact result (IB3I) and I3 
reads in dimensional regularization 

/> = _ 27 1t «r (1 -y)^ ((4 _ <0 (C1) 

For the other integrals, we first carry out the integral, then Taylor expand the result in 
k 2 /o; 2 and finally perform the sums with the zeta function, thus obtaining 

T 2 T(2-d/2)^(^T) d ^ T(d/2-l)T(l + l) (k ^' 

1 "" 12 2tt 2 T ^T(d/2-l-l)T(2l + 2y { ) \2nT t 

(C2) 

k 2 r(2 - d/2)^ ( ^rf ^ r(d/2-i)r(z + i) /* \ a , . 

/2 " 8^ ^ r(d/2 - 1 - l)T(2l + 2) C(2/ + 4 " d) \WTJ ' (C3) 

_ k 4 r( 2 - d/2)^ yzrf f. v_m - i)r(/ + 1) / a, \ 2 < 

/4 " 32^ fa r(d/2 - 1 - Z)r(2/ + 2) C(2/ + 6 " d) V 2^T J ' (C4) 

In the fermionic, sector we have 

rp2 

lo = -24, (C5) 

and we can derive the expansions for I\ and -Z2 following the same procedure used for the 
bosonic integrals, but ending up with the Hurwitz zeta function as a result of the odd 
frequency sums. Thus we have 

- _ r(2 - d/2)/z* (V^T) d ^ T(d/2-l)T(l + l) ( k y 

h ~ 2^r fa r(d/2-i-/)r(2/ + 2) C(2/ + 2 - rf ' 1/2) \wr) ,(C6) 

~ k*r(2-d/2)^(^T) d ~ r(d/2-l)r(Z + l) / fc \ 2 ' 

/2 - 8^T^ JL r(rf/2-i-/)r(2/ + 2) C(2/ + 4 " rf ' 1/2) V^r J • 
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(C7) 
Plugging these expressions in Eqs. (E]) and (I14p we obtain the high-temperature expansion 

(HHD- 

We consider now the low-temperature expansion. The vacuum part gives the order k 2 
term in the expansion, whereas, for the matter part, the condition |k| 3> T translates in 
Eq. (fTTj) into |k| 3> |p|, since the internal momentum |p| is of order T. Expanding this 
expression in |p|/|k| <C 1 yields 

< S (|k| » T) mat = ~9 2 C A ^ + g 2 T 2 H^\ . (C8) 

The singular term (oc |k| 3 /T) and the subtracted zero-mode part also contribute in this 
region. The sum of Eq. (1C8P with the vacuum, subtracted zero-mode and singular parts 
yields Eq. (!2Tj) . For what concerns the static modes, the only scales are |k| and mp, thus the 
condition |k| ^> T becomes |k| ^> irtr, and we end up with Eq. (125]) . Finally, the fermionic 
contribution is suppressed in this region, i.e. the first nonzero term in the expansion of Eq. 
(USD is of order g 2 T 4 /k 2 . 

Appendix D: Non-static two-loop sum-integrals 

We set on the evaluation of the two-loop sum- integrals defined by Eq. ( 1381) . Jo does not 
contribute in dimensional regularization because the integral over k has no scale. Ji can be 
rewritten as 

Jl V J {27T)d J p k , p2q2 /x J {2n)d ^] p q 2 A* J {27r)d J p kVg2 " I ) 

The first term vanishes in dimensional regularization, whereas the second one yields 15 



(27r) d J k 4 pV 8(4 
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Jo can be read from 
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J2 = /ie/ (2^XkV? = (^FV ^ + 1 V~2 + T-^*' (D3) 



15 A convenient way to proceed is by performing first the momentum integrations, by means of two Feynman 
parameters, and then the frequencies sum, which gives £(0) = —1/2. 

44 



J% vanishes in dimensional regularization because the k integral has no scale and finally J4 
yields 






An 



(47T) 



We consider now the fermionic integrals. Jo vanishes because it has a scaleless k integra- 
tion, whereas J± can be computed along the lines of its bosonic counterpart, performing the 
sum over odd frequencies by means of the the generalized (Hurwitz) zeta function, 



J 1 = jj, 



2( 



d d k 

(2tt) c 



r£ 



,2c 



/'" 



d d p to. 



~,2 

n 



TC(0,l/2) 



(2vr) d k 4 pV 4(4vr) 



0. 



Jo can be read from 



50J, 



J 2 = At 



2r 



(2vr) t 



^E" 2 ' 



(i d p 



T 



(2vr) d k 2 pV (4tt) 



In 2. 



(D5) 



(D6) 



Appendix E: Static-modes contribution to the Polyakov loop 



In this appendix, we evaluate the 6-dimensional two-loop integral entering Eq. (T43l) . We 
will perform the calculation modifying the magnetostatic propagator in Eq. (j3]) into 



^(^ 



;i-0^1'^o : 
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(El) 



k 2 / "" k 2 + m*, V tJ k- 

where m m may be interpreted as a small magnetic mass to be put to zero at the end of 
the calculation. The magnetic mass modifies the static gluon self-energy expression with 
resummed gluon propagators from Eq. ( )23i) to 

((k + q)-p) 2 ' 
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where q = k — p. 16 

In Eq. (|43l) . the integral over the first term in Eq. (1E2j) . i.e. the tadpole contribution, 
gives 
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We use here a different parameterization of the integrand with respect to Eq. ([23 
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For the second term, we start by considering the term proportional to (k + q) 2 . We rewrite 



(k + q) 2 = 2(k 2 + m 2 D ) + 2(q 2 + m 2 D ) - (p 2 + m 2 m ) + (m 2 m - 4m 2 ,) 



(E4) 



and consider the contributions given by each of the four terms in brackets. The first one 
gives 
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the second one gives 
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the third one gives 
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and the last one 
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Finally, we consider the term proportional to ((k + q) • p) 2 /p 2 in Eq. (1E2I) . We rewrite the 
numerator as 

(1 - ((k+ p q 2 } ' P)2 = ^[(k 2 + m 2 D ) 2 + (q 2 + ml) 2 - 2(k 2 + m^q 2 + m 2 D )]. (E8) 



The first term gives 
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the third term is —2 times this one and the second term gives 
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The static contribution is thus 
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2 2mo(2mD + m m ) 

The final result is independent of the gauge parameter £. The expression is well behaved for 

m m ->■ and yields Eq. OH]) . 



Appendix F: The Polyakov loop in Feynman gauge 



In this section, we sketch the computation of the vacuum expectation value of the 
Polyakov loop in Feynman gauge. We restrict ourselves to the fundamental representa- 
tion (L = Lp). Since the fermionic contribution, evaluated in Sec. IHIl is to that order 
gauge-invariant, we do not need to compute it here again. 

The perturbative expansion of the Polyakov line through the Baker-Campbell-Hausdorff 
formula is, following 24j and up to order g , 
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and A°(r,x) = A°(r). We recall that 
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where in Feynman gauge the free temporal-gluon propagator is 
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a) b) c) 

FIG. 8: Diagrams contributing to the Polyakov loop up to order g 4 in Feynman gauge. The blob 
stands for the one-loop gluon self energy, the solid line for the Polyakov line and the dots at its 
beginning/end represent the points (0,x) and (1/T, x), which are compactified by the periodic 
boundary conditions. When integrating over loop momenta of order mp, the dashed lines stand 
for resummed temporal propagators, elsewhere for free ones. 



We can now start working on the different terms in Eq. (jFip . The first one gives 
i A.*?* 2 \ ^9 2 C F 2e f d d k 

Following the same approach as in Sec. IIIIl at order g A , the relevant diagrams contributing 
to (1F4[) are shown in Fig. [8] a) and b). At leading order, the Debye mass is gauge invariant, 
whereas the full one- loop gluon self-energy is not. It is convenient to separate non-static 



from static modes. The former yield 49 1 

< s (0, k) = -2g 2 C A {^h ~ (d - l)h + 7 2 ) , (F5) 

where the master integrals Ij are those denned in Eq. ([7]), hence Eq. (IF5I) equals the first 
three terms of the static-gauge expression §6§. The static mode contribution to the self 
energy is common to all gauges that share the same static propagator as the static gauge 
and the Feynman gauge do. Therefore, the static part of the self energy in Feynman gauge 
is just Eq. (|23l) with £ = 1. We then have, separating the contributions coming from the 
scale T from those coming from the scale mp, 

i ng?(|k|~T) n^ki) 



/i 2e / 7 __D 00 (o,k) = A i 2 * 



(2ir) d uuv ' ' ^ J (2tt) 



+ ..., (F6) 



_k 2 + m 2 D k 4 {k 2 + m 2 D ) 2 _ 

where the dots stand for higher orders in the perturbative expansion. We have omitted the 
non-static contribution at the scale m D (cf. Eq. fj40|) ) since it can be shown that also in 

48 



Feynman gauge n^ s (|k| ~ win) — vrvjj = 0(g 2 k 2 ), leading to a higher-order contribution, 
whereas the contribution of the static modes at the scale T leads to a scaleless integral. 
Plugging Eq. (IF6[) into Eq. (]F4j) and using the results of appendices [D] and [E] we obtain 
most of the final, order g 4 , result, except for the contribution of J 4 in Eq. ( 157)) . 

We then consider the other terms in the Baker-Campbell-Hausdorff expansion, starting 
from Hf: 

— /Tr— if 2 \ = F A g 4 / dn / dr 2 / dr 3 / rfr 4 [D o(r 2 - r 3 )Doo(ri - r 4 ) 

--Doo(t 2 - r 4 )-Doo(ri - r 3 )] 

= - 2C(0) 1(£f + --- = ^^ + ---' (F7) 

where we have used free propagators and the dots stand for higher orders. This result 
corresponds exactly to the contribution of J 4 in the static gauge. The contribution can be 
;raced back to diagram c) in Fig. |8] and it corresponds to the term called £ 4 in Eq. (4) of 



2l|. 



We now need to show that the sum of the remaining terms yields zero at order g 4 . 
(Ti2gHoHi) vanishes because it involves a three temporal-gluon vertex. (TY2g 2 HoH2) is a 



more complicated object, however one can show that, working with free propagators 24], 

^(Trg 4 H H 2 ) = + 0{g\g 4 x (m D /T)). (F8) 

The Hq term vanishes, again due to the three temporal-gluon vertex and the H^Hi term 
can be easily shown to be zero after performing the colour trace. The Hq term gives 

^<TY 9 X> = £ (3CJ - ^±) i (> / 1^„„(0, k)) 2 , (F9) 

which is at least of order g 4 x (ttid/T) 2 . This, finally, shows that the Feynman-gauge 
computation of the Polyakov loop agrees with the static-gauge computation that led to 
Eq. 015]). 

Appendix G: Octet contributions 



In this appendix, we want to prove that, up to order g%rT)°, <f )NS = 6°P™\ v .„ Vo , 

5 0(r*)S = _ 6 0(r 2 )S^ ^ ^£ pNRQCD = _^ pNR QCD ) where ^ ]eft _ &nd right . hand sides of 

the equalities encode non-zero modes, zero-modes and higher-multipole one-loop corrections 
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to the pNRQCD octet and singlet propagators respectively induced by interaction vertices 
of the type SV 1 ...r in c\...G\ i _ 1 -E ,ln O + Hermitian conjugate or O^r 11 ...r tn d il ...d in ^ 1 E' ln O + 
charge conjugate. 





=B B= 



a) b) 



/ 






\ 


/ 


\ 


/ 


\ 


A 


\ 


1 


^ 



c) d) 

FIG. 9: The pNRQCD Feynman diagrams giving the leading-order correction to 5 . The single 
continuous line stands for a singlet propagator, the double line for an octet propagator, the circle 
with a cross for the chromoelectric dipole vertex proportional to Va in the Lagrangian (|74p . the 
square with a cross for the chromoelectric dipole vertex proportional to Vb in the Lagrangian (174|) . 
the circle with a dot for the chromoelectric dipole vertex proportional to Vc in the Lagrangian (|74|) . 
the curly line for a chromoelectric correlator and the dashed line for a temporal-gluon propagator. 

The general argument goes as follows. Let's first consider contributions coming from the 
non-zero modes of the loop integral, Fig. [7] providing the leading-order contribution to the 
singlet propagator and diagram a) in Fig. [9] providing the leading-order contribution to the 
octet propagator. As the leading-order example shows, there is a one to one correspondence 
between diagrams in the singlet and in the octet channel, to each singlet diagram corresponds 
an octet diagram whose contribution is equal to the singlet diagram contribution with V s 
replaced by V and viceversa. We note that, since at order g A these contributions are linear 
in V — V s , they are at that order one the opposite of the other. 

Let's now consider contributions coming from the zero modes of the loop integral. In 
order to see how things work, we consider, first, the order r 2 contribution. In the singlet 
channel, only one diagram, Fig. [7J contributes; that contribution has been written in Eq. 
(187}) and evaluated in Eq. ( |88l) . In the octet channel, four diagrams contribute, which are 
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shown in Fig. [9j Diagram a) gives the same contribution as the singlet channel: 

<^ S = -9 2 ^ r ^j^AE^U ab E^(0,k)\\^T + O(g%rT)). (Gl) 

Diagram b) is like diagram a) with the colour factor 1/(2N) replaced by d abc d abc /[4:(N 2 — 1)]: 

C = V^^/^(^^^ 6 )(0,k)| |k| . r + O(, 6 (rT)). (G2) 

Finally, diagrams c) and d) are like diagram a) with the colour factor 1/(2 AT) replaced by 

f abc f abc /[8(N 2 - 1)]: 

Cr 3 = 9 2 j ^/^(E-U^^OM^T + OiArT)), (G3) 

where the positive sign comes from moving a derivative acting on the chromoelectric field 
in one vertex to the temporal gluon in the other one (see also footnote [T2|) . Summing Eqs. 
(IGlj) - (IG3j) we obtain the opposite of the singlet contribution in Eq. (jS] 



This argument may be easily generalized to any order in the multipole expansion. Let's 
consider diagrams contributing to order 2n in the multipole expansion. The singlet contri- 
bution is proportional to 

;OWs 2„iy 1 1 (GA) 



2N ^ (21 + 1)! (2n - (21 + 1))! 



Again there are three classes of octet contributions that correspond to the three classes 
discussed at order r 2 . Except for the first class, each one has a different colour factor with 
respect to the singlet contribution, but for the rest they are equal: 

n-1 

°o,t °c 2Nj^(2£+l)\(2n-(2£ + l))V [ ' 



5?t oc -r 2n 



T 2 a n-l 



N 1 -4 



1 1 

1^ (Off a. A\\ (Or, - (OP A- TW' ^ ' 



AN j^ (2£+l)\(2n-(2£ + l))V 



6^ d)S oc r 2n 






4^(2£)!(2n-2£)! 



where the positive sign in the last expression comes from moving an odd number of deriva- 

n 1 1 
fives acting on the field in one vertex to the field in the other one. Since \^ , , . — 



, q (2£)\ (2n - 2£)\ 
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n-1 



E-r~ n Til i— n ttt, the sum of all octet contributions is just the opposite of the 



singlet contribution. 
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